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Binaries discovered by the SPY project. III. HE 2209-1444: a 
massive, short period double degenerate 1 
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Abstract. In the course of our search for double degenerate (DD) binaries as potential progenitors of type la 
supernovae with the UVES spectrograph at the ESO VLT (ESO SN la Progenitor surveY - SPY) we discovered 
HE 2209— 1444 to be a double- lined system consisting of two DA white dwarfs. From the analysis of the radial 
velocity curve we determined the period of the system to be P = 6 h 38 m 47 s . The semi-amplitudes for both indivi- 
dual components are 109 kms -1 each. A model atmosphere analysis enabled us to derive individual temperatures 
for both components (8490 K and 7140 K, resp.) and masses of O.58M0 for each component. The total mass 
of the system is 1.15 ± O.O7M0. The system will loose angular momentum due to gravitational wave radiation 
and therefore will merge within 5Gyrs - less than a Hubble time. HE 2209— 1444 is the second massive, short 
period double degenerate detected by SPY. Its total mass is about 20% below the Chandrasekhar mass limit and 
therefore it does not qualify as a potentional SN la progenitor. However, together with our previous detections it 
supports the view that Chandrasekhar mass systems do exist. 
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1. Introduction 

Supernovae of type la (SN la) play an outstanding role 
in the study of cosmic evolution. In particular, they are 
regarded as one of the best standard candles for the de- 
termination of the cosmological parameters Hq, f2, and A 
(e.g. Riess et al. 119981 Leibundgut I2001|l . However, the 
nature of their progenitors remains a mystery (e.g. Livio 

There is general consensus that the SN event is due 
to the thermonuclear explosion of a white dwarf likely 
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when the Chandrasekhar limit of w 1.4 Mq is reached. But 
the exact nature of the progenitor system is still unclear. 
While it must be a binary with matter being transferred 
to the white dwarf from a companion, two main scenar- 
ios exists. According to the so-called double degenerate 
(DD) scenario (Iben and Tutukov ll984fl . the mass donat- 
ing companion is also a white dwarf. In the so-called singe 
degenerate (SD) scenario (Whelan and Iben 1973), the 
mass donating component is a red giant /subgiant. 

A DD has to fulfill two criteria to be considered as 
a potential SN la progenitor. The total mass of the sys- 
tem has to exceed the Chandrasekhar limit, and the two 
components have to merge in less than a Hubble time 
due to the loss of angular momentum via gravitational 
wave radiation. During the last decades, several system- 
atic radial velocity (RV) searches have been undertaken 
(see Marsh [2000 for a review). By now, combining all 
surveys nearly 180 white dwarfs have been checked for 
RV variations yielding a sample of 18 DDs with periods 
of less than 6.3 days (Marsh I2U0"U1 Maxted et al. l2000afl . 
Among these 18 systems, only 6 double- lined systems have 
been found (in three of them, the companion is barely de- 
tectable). However, none of the known 18 systems seems to 
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be massive enough to qualify as a SN la precursor. Since 
theoretical simulations suggest only a few percent of all 
close DDs to be potential SN la progenitors (Iben et al. 
1997, Nelemans et al. I2()()l|l this result is not surprising 
and a larger sample size is needed to find a potential SN 
la progenitor. 

Therefore -to perform a definite test of the DD 
scenario -we have embarked on a large spectroscopic 
survey of ks 1000 white dwarfs using the UV- 
VisualEchelle Spectrograph (UVES) at the ESO VLT 
UT2 (Kueyen), searching for RV variations of white dwarfs 
and pre white dwarfs (ESO SN la Progenitor surveY 
- SPY). The ongoing SPY project already provided a 
wealth of new RV variable DDs (Napiwotzki et al. |2001a 
and 2003). In particular, more than a dozen double- lined 
systems have been found by SPY. These are of outstand- 
ing importance because the radial velocity curves can 
be solved for all system parameters if the gravitational 
redshifts of the components are unequal (see Sect. UJ. 
An analysis of the subdwarf B plus white dwarf system 
HE 1047—0436 was presented in the first paper of this se- 
ries (Napiwotzki et al. [2001b) . In paper II (Napiwotzki 
et al. I2002[) we have presented the double-lined system 
HE 1414-0848, which has a system mass of 1.26 A/©, but 
will merge only in two Hubble times. 

In this paper we report on the results of our 
follow-up spectroscopy of the new double-lined system 
HE 2209-1444. We will show that this system is another 
massive DD, but with an orbital period shorter than the 
HE 1414—0848 system. After a description of the obser- 
vations and data reduction in Sect. |21 we deal with the 
determination of the system's radial velocity curves and 
orbital parameters in Sect. [21 Section 0] gives a model at- 
mosphere analysis of the spectra. We compare the results 
to predictions of population synthesis models in Sect. [3] 
and summarize in a final Sect. 

2. Observations and data analysis 

HE 2209-1444 (a 20 oo = 22 h 12 m 18.1 s , <5 2000 = -14°29'48", 
Bpg = 15.3 mag ) was discovered by the Hamburg ESO sur- 
vey (HES; Wisotzki et al. 12213 Christlieb et al. lMlTjl as a 
potential cool white dwarf and, therefore, was included in 
our SPY project. From the first survey spectrum taken in 
December 2000 we found HE 2209-1444 to be a double- 
lined binary system consisting of two DA white dwarfs. 
The H Q line cores of both components were separated by 
4.5 A in the discovery spectrum, corresponding to a RV 
difference of 200km/s. 

Therefore we observed the system during our follow- 
up observations for the SPY project at the Calar Alto 
Observatory, Spain. From 2001 July 6 th to July 15 th a to- 
tal of 15 medium resolution spectra were obtained with the 
TWIN spectrograph at the 3.5 m telescope. For the blue 
arm we used grating No. 5 while No. 6 was used in the red. 
Both gratings have 1200 lines/mm, giving a dispersion of 
0.54A/pix. The wavelength coverage ranges from 3900 A 
to 5000 A at a resolution of 1.3 A in the blue part of the 
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Fig. 1. Spectra of HE 2209-1444 taken at the 3.5 m tele- 
scope of the Calar Alto Observatory during one night. The 
observed spectra are as well plotted as the fits. The JD 
is heliocentric corrected and computed for the mid of the 
exposure. 



spectrum. In the red part, the wavelength scale ranges 
from 6000 A to 7000 A at a resolution of 1.2 A. The data 
were reduced with a set of standard routines for long slit 
spectra from the ESO MIDAS package. A two dimensional 
bias correction was done for the object data as well as for 
the dome flats. After flatfield correction, the spectra were 
extracted from the two dimensional images with respect 
to the sky background. Finally we performed the wave- 
length calibration using the internal ThAr lamp of the 
spectrograph. 

Figure ^ shows a sequence of six Hq, spectra taken at 
Calar Alto Observatory during two hours. The observed 
spectra are plotted as well as the fits (cf. Sect. |3J). The 
rapid change of the spectral appearance due to the orbital 
motion is obvious. 
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Our TWIN spectra are supplemented by spectra taken 
with the ESO VLT and the Isaac Newton Telescope. The 
UVES data we got from the ESO VLT UT2 were taken in 
December 2000 and October 2001. Although the data pro- 
vided by the UVES reduction pipeline were quite useful 
for candidate selection, in many cases the quality of the 
pipeline data was insufficient for a quantitative analysis 
of Balmer line profiles. Therefore, the raw data were re- 
duced with a semi-automatic reduction package developed 
in Bamberg (details will be published in a forthcoming pa- 
per). 

The last set of data used for our analysis of 
HE 2209-1444 was taken on October 28 th to October 30 th 
using the 2.5 m Isaac Newton Telescope (INT) on La 
Palma. Eight spectra were taken with the Intermediate 
Dispersion Spectrograph (IDS), using the 500mm cam- 
era and R1200R grating, giving a dispersion of 0.37A/pix 
from 6300 to 6700 A at a resolution of 1 A. The spectra 
were reduced (using the overscan bias level and a tung- 
sten flat field) with standard IRAF 1 tasks, using optimal 
extraction fHorne H986|l . The wavelength calibration was 
done using CuAr and CuNe lamp taken at the same posi- 
tion of the telescope and extracted at the same place on 
the CCD as the object spectra, typically giving an rms 
scatter less than 0.01 A. 



3. Radial velocity curve and orbital parameters 

RVs were determined from the TWIN and UVES spec- 
tra from a fit of a set of mathematical functions to the 
observed line profiles within the MIDAS context. We con- 
centrated on a 100 A wide range of the spectra centered 
on H a . Two Gaussians were used to fit both line cores, a 
Lorentzian to model the line wings and a linear function to 
reproduce the overall spectral trend. Radial velocities of 
both components were calculated from the central wave- 
lengths of the fitted Gaussians. Afterwards the measured 
RVs were corrected to heliocentric values. The INT spec- 
tra were analyzed using the MOLLY software package, 
fitting three Gaussians to the central part of the H a line. 

Both components are very similar, but the left one in 
Figure^is slightly deeper and broader. This component is 
called component A further on. However, sometimes it was 
difficult to do an unambiguous identification of the com- 
ponents especially near conjunction phases and in some 
spectra with relatively low S/N ratio and of lower reso- 
lution (from Calar Alto and La Palma). Therefore, we 
computed the separation of the components first and fit- 
ted sine curves for a range of periods to the measured 
RV differences and determined the x 2 value for each pe- 
riod. All times were heliocentric corrected and calculated 
for the center of the exposure. Doing this we produced a 
" difference power spectrum" indicating the quality of our 
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RV fit as a function of period (like the one presented in 
Fig. 0). This way we got a first estimate of the period of 
the system and were subsequently able to identify the in- 
dividual components unambiguously for all of our spectra. 
This enabled us to compute for each individual component 
a "power spectrum". Both power spectra showed an out- 
standing peak at 6 h 38 m 46.2 s (0.276924 d ). An inspection 
of the phased RV curves created using periods correspond- 
ing to other peaks in the power spectra allowed us to rule 
out aliases, because of one or more strongly deviating RV 
measurements. This way, the period corresponding to the 
main peak of the power spectrum remained as an unam- 
biguous solution for component A as well as for compo- 
nent B - exactly as it should be for a binary system (see 
Fig. • Finally, adding the \ 2 values of the two individual 
power spectra we produced a " combined power spectrum" 
for the whole system (Fig. ■ 

In a final step we determined accurate orbital param- 
eters from a simultaneous fit of all spectra with the pro- 
gram FITSB2 (see Sect.0]and Napiwotzki et al.[2003) in a 
fashion similar to the method described in Maxted et al. 
(2002). Two model profiles, one for each star, consisting 
of a combination of two Gaussians for the line core and a 
Lorentzian for the line wings were fitted. The position of 
each profile is determined from 



RV A /b = 7a/b ± K A/B sin 



T-Tp 
P 



Free parameters are the mean velocities (including gravi- 
tational redshift) of each component 7a/Bj the (projected) 
orbital velocities -f^A/B, the zeropoint To, and the period 
P, plus the parameters defining the line profiles. 

The final ephemeris of the system for the time To de- 
fined as the conjunction time at which star A moves from 
the blue side of the RV curve to the red one (i.e. star A is 
closest to the observer) is 

Hel.JD(To) = (2, 452, 097.0349 ± 0.0021) + 

(0.276928 ± 0.000006) x E. (1) 

Also from the RV curves (see Fig. |3J| of the FITSB2 
analysis we found the semi-amplitude of component A 
to be K A = 108.7 ± 7.5km/s (7A = -14.4 ± 4.5 km/s) 
whereas for component B the semi-amplitude is Kb = 
109.1 ±10. 75 km/s ( 7B = -18.7 ±6.8 km/s). The mass ra- 
tio can be computed from the ratio of the semi-amplitudes 
to be 



Mb 
M a 



Ka 



1.00 ±0.13 



(2) 



i.e the masses of both white dwarfs are identical within 
the error limits. For white dwarfs the gravitational red- 
shift can easily be measured from high resolution spec- 
tra (e.g. Reid ll99rj|l . As described in paper II for the case 
of HE 1414—0848, in a double degenerate system the dif- 
ference of gravitational redshifts can be measured from 
7A — 7b ■ If this difference is non zero we obtain an addi- 
tional constraint equation, which allows to solve for the in- 
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Fig. 2. Combined power spectrum of the system 
HE 2209-1444. The x 2 values obtained from the individ- 
ual power spectra of the individual components are coad- 
ded. 
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Fig. 3. Radial velocity curve of the system 
HE 2209-1444. The measured RVs for component A 
(filled triangles) and B (open triangles) are plotted 
against the orbital phase. 

dividual masses. However, since the masses of both compo- 
nents of HE 2209—1444 are very similar, the redshift differ- 
ence is zero within error limits. Hence we cannot solve for 
the individual masses from the radial velocity curve alone. 
The required additional information, however, will be ob- 
tained from a quantitative spectral analysis described in 
the next section. 

4. Spectroscopic analysis 

Since in the case of HE 2209— 1444 mass estimates from 
the RV curves alone are impossible, we had to rely on a 
model atmosphere analysis. Because this system is double- 
lined the spectra are a superposition of both individual 



white dwarf spectra. This problem is sometimes circum- 
vented by assuming that the results of a simple single- lined 
model fit is representative for the average parameters. In 
paper II we estimated the mean gravity and individual 
temperatures of the components of HE 1414—0848 from 
the analysis of individual spectra taken close to conjunc- 
tion and quadrature phases. However, this approach is of 
limited accuracy, only. 

A direct approach would be to disentangle the ob- 
served spectra by deconvolution techniques into the spec- 
tra of the individual components. Then we could ana- 
lyze the spectra by fitting synthetic spectra developed for 
single-lined white dwarfs to the individual line profiles. 
Such a disentangeling procedure has been developed by 
Simon and Sturm l|1994|l and was sucessfully applied to 
main sequence double-lined binaries. However, it has not 
been tested for white dwarfs, for which the wavelength 
shifts caused by orbital motions are much smaller than the 
line widths of the broad Balmer lines. Therefore we choose 
a different approach for our analysis of the HE 2209— 1444 
system. We used the program FITSB2 (Napiwotzki et al. 
2003 ) , which performs a spectral analysis of both compo- 
nents of double-lined systems. It is based on a % 2 mini- 
mization technique using a simplex algorithm. The fit is 
performed on all available spectra covering different spec- 
tral phases simultaneously. For HE 1414—0848 (paper II) 
we made use of spectra taken at conjunction or quadrature 
only. But now, all available spectral information is com- 
bined into the parameter determination procedure. An ap- 
plication of the program fitsb2 to HE 1414—0848 yielded 
encouraging results (Napiwotzki et al. 120030 . A detailed 
description of the program and an evaluation of its per- 
formance will follow in a subsequent paper (Napiwotzki, 
in prep.). 

A large grid of synthetic spectra for DA white dwarfs 
computed from LTE model atmospheres with a code de- 
scribed in Finley et al. 1)1997(1 was used for the analysis. 
A simultaneous fit of the Balmer lines H^ to Hg (UVES 
spectra) or Hp to H$ (TWIN spectra) was performed. 
We did not include the INT spectra for the model at- 
mosphere analysis, because these spectra cover only the 
H Q range. For details refer to Koester et al. I|2001l and ref- 
erences therein). The model spectra were convolved with 
Gaussians with FWHMs corresponding to the resolution 
of the observed spectra. 

The total number of fit parameters (stellar and or- 
bital) is high. Therefore we fixed as many parameters as 
possible before performing the model atmosphere analy- 
sis. We have kept the radial velocities of the individual 
components fixed according to the radial velocity curve 
presented in Sect. [21 Since the mass ratio is already accu- 
rately determined from the radial velocity curve we fixed 
the gravity ratio. The remaining fit parameters are the ef- 
fective temperatures of both components and the gravity 
of the primary. The gravity of the secondary is adjusted ac- 
cording to the primary value during the fitting procedure. 
The surface gravities also determine the relative weight 
of the two model spectra from the radius, obtained from 
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Table 1. Fitted heliocentric radial velocities. The heliocentric corrected JDs are given for the center of the exposures. 
The phases <f> are computed from the ephemeris given in Eq. ^ At $ = 0,982 and $ = 0,040 the individual RVs of 
both components can not be measured accurately, because the spectra were taken close to conjunction. However, these 
spectra were used to determine the orbital parameters using the procedure described in Sect.OH 



Hel.JD 


$ 


Heliocentric RV [km/s] 


Telescope 


Instrument 


-2,400,000 




A 


B 






51885.5488 


0.314 


85 ± 3 


-119 ± 4 


ESO UT2 


UVES 


52097.5404 


0.826 


— 105 ± 6 


55 ± 9 


CA 3.5 m 


TWIN 


52097.6125 


0.086 


67 ± 8 


—36 ± 10 


CA 3.5 m 


TWIN 


52101.5702 


0.377 


63 ± 5 


—91 ± 8 


CA 3.5 m 


TWIN 


52101.6309 


0.596 


-77 ± 4 


59 ± 6 


CA 3.5 m 


TWIN 


52102.5410 


0.883 


-78 ± 8 


68 ± 13 


CA 3.5m 


TWIN 


52102.5685 


0.982 






CA 3.5 m 


TWIN 


52102.5845 


0.040 






CA 3.5 m 


TWIN 


52103.5514 


0.531 


-51 ± 10 


27 ± 12 


CA 3.5 m 


TWIN 


52103.5678 


0.591 


—89 ± 12 


41 ± 14 


CA 3.5 m 


TWIN 


52103.5850 


0.653 


— 100 ± 6 


55 ± 8 


CA 3.5 m 


TWIN 


52103.6014 


0.712 


—124 ± 5 


88 ± 7 


CA 3.5 m 


TWIN 


52103.6182 


0.773 


—120 ± 4 


83 ± 7 


CA 3.5 m 


TWIN 


52103.6344 


0.831 


-97 ± 4 


77 ± 8 


CA 3.5 m 


TWIN 


52103.6504 


0.889 


—82 ± 5 


45 ± 9 


CA 3.5 m 


TWIN 


52103.6666 


0.948 


-41 ± 8 


37 ± 11 


CA 3.5 m 


TWIN 


52104.6750 


0.589 


-85 ± 8 


73 ± 12 


CA 3.5 m 


TWIN 


52105.6591 


0.143 


63 ± 5 


-103 ± 7 


CA 3.5 m 


TWIN 


52106.5183 


0.245 


99 ± 6 


-108 ± 7 


CA 3.5 m 


TWIN 


52172.6372 


0.004 


-14 ± 2 


-14 ± 2 


ESO UT2 


UVES 


52194.5343 


0.075 


35 ± 3 


-70 ± 5 


ESO UT2 


UVES 


52194.6249 


0.402 


47 ± 4 


-84 ± 6 


ESO UT2 


UVES 


52211.3592 


0.831 


-121 ± 5 


36 ± 7 


INT 


IDS 


52211.3807 


0.908 


-76 ± 5 


26 ± 8 


INT 


IDS 


52212.3189 


0.296 


72 ± 5 


-112 ± 7 


INT 


IDS 


52212.3440 


0.387 


31 ± 4 


-100 ± 6 


INT 


IDS 


52212.3921 


0.561 


-85 ± 4 


18 ± 7 


INT 


IDS 


52212.4138 


0.639 


-114 ± 4 


32 ± 6 


INT 


IDS 


52212.4822 


0.886 


-106 ± 7 


26 ± 11 


INT 


IDS 


52213.4533 


0.393 


32 ± 5 


-88 ± 8 


INT 


IDS 



the mass-radius relation of Benvenuto & Althaus (1999). 
The flux ratio in the V-band is calculated from the actual 
parameters and the model fluxes are scaled accordingly. 
The individual contributions are updated consistently as 
part of the iteration procedure. 

Strong NLTE cores are present in H a and H^, which 
cannot be reproduced by our LTE model spectra. In prin- 
ciple, these lines are important for the temperature de- 
termination, especially of the cooler and fainter B com- 
ponent. However, due to the large NLTE effects we had 
to exclude H Q completely and the core of H^ (± 4A). 
The final fit results are summarized in Table |21 and a 
sample fit is shown in Fig. The fit of the line cores 
is not perfect, even after excluding H a . The parameter 
which is most sensitive to details of the fitting procedure 
is, as may be expected, the temperature of the fainter 
component B. If, e.g., we include the wings of the H Q line 
we increase T e g of component B by almost 500 K, while 
the temperature of A and the gravity are only modified 
within the error limits. Thus, while the temperature of B 
is less well determined than that of A, the basic proper- 



ties of the HE 2209— 1444 system are reliable. Both com- 
ponents have identical masses (Ma = Mb = 0.58 M Q ) but 
different temperatures (8490 K vs. 7140 K). White dwarf 
masses were computed from the mass - radius relations of 
Benvenuto & Althaus lfTU99) with "thick hydrogen en- 
velopes" (Mh/Mwd = 10" 4 ). This result does not de- 
pend much on the choice of a particular model computa- 
tion. From the thick envelope cooling sequences of Wood 
(19951) we derived virtually identical results: 0.58 M and 
0.57M Q , respectively. However, since the HE 2209-1444 
system is the result of a common envelope evolution it is 
not clear which envelope layer is the correct one. Using the 
"thin" envelope models (Mh = 0) of Benvenuto & Althaus 
(1999} as the other extreme would yield slightly lower 
masses (0.57M© and 0.56 M©, resp.). The resulting sum 
of masses of the HE 2209-1444 system is 1.15 ± 0.07M©. 
The error limit for the sum of masses is smaller than ex- 
pected from a simple combination of the individual errors, 
because the mass errors of the individual components are 
partially anti-correlated. 
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Fig. 4. Sample Model atmosphere fits of HE 2209— 1444 phase spectra taken at different phases: The discovery 
spectrum ($ = 0.31) taken with the UVES spectrograph and two TWIN spectra ($ = 0.60 and $ = 0.71). 



Table 2. Results of the model atmosphere analysis with 
FITSB2. The error limits include the uncertainties of or- 
bital parameter determination. 



comp. 




logs 


M/Mq 


age/Gyrs 


A 


8490 ± 80 


7.97 ±0.05 


0.58 ±0.03 


0.9 


B 


7140 ±110 


7.97 ±0.13 


0.58 ±0.08 


1.4 



5. Evolutionary status 

From the temperatures derived above we can estimate 
the ages of both components using cooling tracks for 
white dwarfs. Our estimate is based on the calculations of 
Benvenuto & Althaus (|1999H for C/O white dwarfs with 
thick hydrogen envelopes. Ages were interpolated from 
these tracks for the derived masses and temperatures of 
the white dwarfs. Resulting ages are 0.9Gyrs for A and 
1.4Gyrs for B, i.e. the cooling ages of both white dwarfs 
differ by 0.5Gyrs or 50%. 



The separation between both white dwarfs is quite 
small, only 1.87i? Q . Thus HE 2209-1444 obviously un- 
derwent phases of strong binary interaction in its history. 
From the size of the orbits and the period the orbital ve- 
locity can be computed to be 171km/s. The comparison 
with the observed RV amplitudes allows us to determine 
the inclination of this system as i = 40° . This system will 
merge within 5 Gyrs due to gravitational wave radiation. 
Figure [S] shows the total masses and periods of the sample 
of known double-lined DD systems. As can be seen, nei- 
ther HE 2209-1444 nor HE 1414-0848 (Paper II) qualify 
as a SN la progenitor. However, both objects found by 
SPY are much closer to the region where such progenitors 
are expected than any DD system known before. 

We can compare the properties of HE 2209—1444 with 
the outcome of computations for the synthesis of the pop- 
ulation of double white dwarfs of Nelemans et al. H2001fl . 
Fig. shows the expected period - mass ratio distribution 
of systems in which one component has a mass above 
0.50 Mq from model A2 of Nelemans et al. lf20"0"T|) . Only 
systems with a luminosity ratio between the components 
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Fig. 5. Mass - period diagram of known double- lined DD 
systems with mass determinations for both components. 
Literature values from Maxted et al. ©99 and 2000c) are 
plotted as open circles while SPY results are shown as 
filled circles with error bars (WD 1349 +144 parameters 
from Karl et al. 2003). The Chandrasekhar limit and the 
region of DDs, which will merge within a Hubble time are 
indicated. 

of less than 5 are plotted. The derived period and mass 
ratio of HE 2209— 1444 is marked in Fig. with its error 
bar. The figure shows that such systems are expected to 
have periods smaller than «10 hours and either similar or 
more massive companions. 

HE 2209-1444 confirms the trend that double white 
dwarfs tend to have roughly equal masses, in contrast 
to population models using a standard common enve- 
lope description (see for a discussion Maxted et al. UOOOc 
Nelemans et al. [2.001 ,) and more in line with the results 
presented in Nelemans et al. 120011 in which the outcome 
of dynamical unstable mass transfer is calculated from the 
angular momentum balance, rather than the energy bal- 
ance. However, it should be noted that a model using stan- 
dard common envelope also produces systems similar to 
HE 2209-1444. 

The SPY project will provide a large, homogeneous 
sample of close double white dwarfs, which will enable 
us to confirm or refuse the alternative common envelope 
formalism. 

6. Summary 

We report the discovery of the close double-lined bi- 
nary HE 2209-1444 by SPY. From follow-up observa- 
tions at three different telescopes we derived the radial 
velocity curve. The analysis yielded an orbital period of 
P = 6 h 38 m 47 s and identical semi-amplitudes of 109km/s 
for both components. Parameters of both individual com- 
ponents were determined from a model atmosphere anal- 
ysis. We made use of a novel method which combines all 
available spectra covering the different orbital phases for 
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Fig. 6. Period- mass ratio distribution for the current 
population of DDs in the Galaxy for systems in which 
at least one white dwarf has a mass above 0.5 M Q from 
model A2 of Nelemans et al. (2001}. Only systems with 
a luminosity ratio between the components of less than 5 
are plotted. Grey shades indicate the expected number of 
systems brighter than V = 15.5 m . The derived period and 
mass ratio of HE 2209—1444 is plotted as a filled dot with 
an error bar. Also plotted are HE 1414- 0848 (paper II) 
and PG 1115+166 (Bergeron et al., HH and Maxted et 

al., 12002). 

a simultaneous fit. The result is a common gravity of 7.97 
and temperatures of 8490 K and 7140 K for components A 
and B, respectively. 

Masses were determined from a subsequent comparison 
with the cooling tracks of Benvenuto & Althaus 1^999). 
The results are masses of O.58M0 for the individual com- 
ponents and a total mass for the HE 2209—1444 system of 
1.15±O.O7M . This makes HE 2209-1444 the second most 
massive DD ever found in a RV survey, surpassed only 
by HE 1414-0848, another detection of the SPY project 
(Napiwotzki et al. I2002f) . However, the seperation of the 
components in HE 2209— 1444 is smaller than in the case 
of HE 1414-0848. Therefore, the merging time of the for- 
mer system is 5 Gyrs whereas in the latter case the system 
will merge within two Hubble times. 

We compared the resulting orbital parameters of 
HE 2209— 1444 with the outcome of theoretical calcula- 
tions for the formation of DDs and find that the properties 
of HE 2209—1444 are in agreement with recent predictions. 
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